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Abstract 

The structural invariant subspaces of the discrete-time singular Hamiltonian system are used in [T) 
to give an analytic nonrecursive expression of all the admissible trajectories. A deeper insight into 
the features of these subspaces, particularly focused on the dimensionality issue, is the object of this 
note. 



I. Main Content 

In JT), the structural invariant subspaces Vi and V2, defined by (17) and (18) respectively, 
are used in (19), that is the analytic nonrecursive expression of the set of the admissible 
solutions of the discrete-time singular Hamiltonian system (10) over the time interval 

0<k<k f -l. 

In this note it will be shown that the dimension of the subspace V2 may be lower than 
n, where n denotes the dimension of the state space of the original system as defined in 
(1), (2). In particular, the possible loss of dimension of V 2 (or, equivalently, the possible loss 
of rank of the matrix V 2 defined by (18)) depends on the properties of the original system 
(1),(2) (under assumptions A.l-AA). 

Let us consider system (1),(2) and perform the similarity transformation T — [Ti T 2 ], 
where imTi = 1Z, the reachable subspace of (A, B). With respect to the new basis, 
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Moreover, the solution P + of the Riccati equation (11), (12), partitioned accordingly, is 
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where P c is the stabilizing solution of the Riccati equation restricted to the sole reachable 
part of the original system: i.e., 

P c = A]P C A C + CjC c - (A T C P C B C + CjD)(D T D + Bj P c B c )-\Bj P C A C + 7J T C C ), 

with 

D T D + BjP c B c > 0. 
The stabilizing feedback is partitioned as 

-K+ =[K C K u ]. 



Similarly, the solution W of the discrete Lyapunov equation has the structure 
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where W c is the solution of the discrete Lyapunov equation restricted to the sole reachable 
part of the original system: i.e., 
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Simple algebraic manipulations, where these partitions are taken into account, yield the 
following structure for the matrix V 2 : 
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where the symbol * denotes a possibly nonzero submatrix. The structure pointed out in the 
partitioned matrix V2 shows that the rank of V 2 may be lower than n. This circumstance 
occurs, for instance, if A u has a zero row, like in the illustrative example considered in the 
following section. 

It is worth noting that, by contrast, the subspace 
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that is used in (20) of flTJ in order to express the sole state and costate trajectories over the 
time interval < k < kf, has dimension n. 

In fact, the corresponding partitioned matrix is 



V 2 



The rank of V 2 is n, since the symmetric positive definite W c , being the solution of the 
restricted Lyapunov equation above, has the same rank of the controllability Gramian of the 
pair (A c + B C K C , B c ), which is completely controllable by construction. 
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II. An Illustrative Example 

This section presents a numerical example where the rank of matrix V 2 is lower than 
the dynamic order n of the original system, while n is the rank of matrix V 2 . The variables 
are displayed in scaled fixed point format with five digits, although computations are made 
in floating point precision. Consider system (1), (2) in (TJ, with 
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By pursuing the procedure illustrated in [1], one gets, in particular 
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III. Conclusions 

In this note, it has been shown that the dimension of the structural invariant subspace V2 
may be lower than the dynamic order n of the original system. The result has been illustrated 
by a numerical example. This is the reason why the only-if part of the proof of Theorem 1 in 
[DO does not rely on a dimensionality count, but on the maximality of the subspace V2 (and of 
the subspace Vi). In fact, maximality follow from Property 2 (and Property 1, respectively), 
according to [2, Section 5.4]. 
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